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Abstract 

The computation of detection probabilities and arrival time dis- 
tributions within Bohmian mechanics in general needs the explicit 
knowledge of a relevant sample of trajectories. Here it is shown how 
for one-dimensional systems and rigid inertial detectors these quanti- 
ties can be computed without calculating any trajectories. An expres- 
sion in terms of the wave function \E' and its spatial derivative d x fy, 
both restricted to the boundary of the detector's spacetime volume, 
is derived for the general case, where the probability current at the 
detector's boundary may vary its sign. 

1 Introduction 

A microscopic object s may trigger a sudden change in the properties of 
a macroscopic system S. Often such events take place during an interval 
of time much shorter than the duration of the interaction between s and 
S. Perhaps the simplest example is the detection time of a quantum parti- 
cle which slowly passes a detector. This phenomenologically well accessible 
quantity is usually termed arrival time. Quantum theory has difficulties in 
identifying events and a fortiori their time of occurrence within its formalism, 
because, according to Schrodinger's equation, the unitary state evolution of 
the closed system containing s and S does not make any sudden jumps. 
Jumps of states are introduced into the standard quantum formalism only 
through state reduction, which is supposed to happen in an open system, 
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when a "measurement" is performed on it from "outside". Yet in this case 
the instant of time of the reduction is not stochastic but rather determined 
by the observer's deliberate choice. Thus state reduction does not seem to 
be the proper notion to understand the stochastic distribution of the time of 
events. 

Attempts to obtain the arrival time distribution through the model of 
continuous observation leads to the well known quantum zeno paradox p. 
Also the proposals for a time operator (see e.g. f|) are still subject to dis- 
cussion and controversy 0111 . 

One strategy to incorporate detection events into quantum theory is by 
means of Bohmian mechanics (BM), which introduces the additional notion 
of particle trajectories into the standard formalism. Within this framework 
Leavens [3] and McKinnon and Leavens jH] derived an expression for the time 
resolved detection probability for one-dimensional (ID) scattering situations. 
Let \I/ be a normalized solution of a time dependent ID Schrodinger equation 
and j its associated probability current density. In |5J|B] it has been argued 
for a Bohmian particle with wave function \l/ that the detection probability 
at position a between time and time T > is given by f Q \j(t, a) \ dt The 
line of argument assumes from the outset that no trajectory passes through 
a more than once during the time interval [0, T] which is guaranteed if a) 
does not change sign. (For an ideal detector the first entry triggers an event, 
e.g. by discharging the device and producing a click. Further on the detector 
is insensitive to additional entries.) 

If, however, multiple crossings do occur, the replacement of j(t,a) by a 
cut off current has been advocated by Daumer et al jS], such that only the 
first traversal of the trajectories should be counted. This means that the time 
intervals with second, third, etc. crossings should be dropped from the inte- 
gral f Q \j(t, a) | dt. Therefore the computation of such detection probabilities 
in general demands the explicit knowledge of the Bohmian trajectories of the 
problem at hand. Yet Bohmian trajectories are the solutions of a nonlinear 
system of ordinary differential equations and therefore difficult to obtain. 

The general Bohmian notion of detection probability associated with quite 
arbitrary space-time regions has been formalized in j^]. As in the earlier 
treatment jH] Bohmian trajectories enter the formula defining the detection 

1 Under certain provisos the quantum optical detection model of reference concurs 
with this expression. 
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probability. Here we show how for ID systems and rigid inertial detectors the 
Bohmian detection probability and the associated arrival time distribution 
can be computed under quite general circumstances without any knowledge 
of the Bohmian trajectories. The relevant formula is contained in proposition 
□ 

In section |2] we give a very short overview of the main ingredients of 
Bohmian mechanics, omitting mathematical detail. Section |3] introduces 
the essential mathematical structures needed for the formulation of one- 
dimensional arrival time problems in the framework of non-relativistic Bohmian 
mechanics. In section 0] then a reformulation of the arrival time distribution 
for one-dimensional detectors occupying spatial intervals [a, b], with the mere 
aid of probability and current density integrals, is given and proved. Finally 
section |5] demonstrates the practical use of the technique given in section 0] 
by means of several numerical examples. Free evolution as well as evolution 
under the influence of external potentials is considered. 

2 Bohmian Mechanics 

BM rests on the insight that with each normalized solution of the time de- 
pendent Schrodinger equation, a fibration of the configuration space-time is 
given. At time t each fibre (Bohmian trajectory) has a unique representative 
in the underlying configuration space, and therefore a dynamical evolution 
of configurations along the fibres follows. The local conservation of con- 
figuration space probability implies that its quantum mechanical evolution 
coincides with the one implied by the transport along the fibres. There- 
fore a causal deterministic interpretation of quantum mechanics in terms of 
movement in configuration space becomes consistent. An individual quan- 
tum system with wave function \l/ within BM is now assumed to realize one 
of the system's trajectories, i.e. at each instant the system is in a point of the 
configuration space. Amending the continuum notions of quantum theory by 
such point structures opens up the possibility to identify unique properties 
and sudden events within the formalism. Consequently one may ask again 
"Does a trajectory enter a certain spacetime region?" or "When is it, that 
the system's trajectory enters a certain spacetime region?" . 

Since the choice among the trajectories is beyond control, only probabilis- 
tic predictions can be made, yet at any time a system has properties without 
the need to invoke state reduction. Whenever the probabilistic predictions 
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of BM can be compared with those of standard quantum mechanics, they 
agree. To us the prime achievement of BM seems to be, that it resolves the 
quantum measurement problem. A concise summary of Bohmian mechanics 
can be found in [TU] . 



3 Arrival time from Bohmian flow 

For the sake of simplicity, the set of Galilean spacetime points is assumed to 
be M. = M 2 . As a positively oriented, global and inertial chart we choose 
id>M ='■ (t,x). The associated tangent frame is denoted by (d t ,d x ). 

3.1 Bohmian velocity vector field and Bohmian flow 

Let the mapping ^ : A4 — > C be C 2 and a solution to the Schrodinger 
equation 



with V any real scalar potential. Then with the aid of the position and 
current densities 

p ■= 

and 

m 

the current vector field j is defined as 

j := pd t +jd x 

on Ai. If p 7^ 0, Vp G M. the corresponding Bohmian velocity vector field 

v := d t + -d x 
P 

is C l . The maximal integral curve of the vector field v through a point p G M. 
is the unique function 7 : I p — > M. with 7(0) = p and 

7 ( s ) = v(j(s)), Vs G I p 

where I p is a non-extendable open real interval. Those integral curves are 
regarded to represent the worldlines of the actual Bohmian particles. We 
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assume that v is complete, which means that I p = IR, for all p G A4. Then 
the mapping 

F : R x M -> A*, ^ 7 P (s) 

is a global flow on .M with F(r, F(s,p)) = F(r + s,p). As t o F(s,p) = 
t o 7p(s) = s + i(p), no worldline begins or ends at a finite time. Moreover 
F(s, •) bijectively maps instantaneous spaces T, T = {p e M / t(p) = r} onto 
instantaneous spaces, namely F(s, S r ) = T, s+r . 



3.2 Conservation of probability 

By inserting j as the first argument into the volume form E = dt A dx, we 
get the current- 1-form 

J := jjE = pdx — jdt. 

Due to the continuity equation dtp + d x j = the current form J is closed: 

dJ = (d t p + d x j)dt A dx = 0. 

For any spacetime region D with piecewise (^-boundary dD, Stoke's theorem 
therefore assures that 

f J= f dJ = Q. (1) 

JdD JD 

In particular for spacetime regions of the form indicated in figure ^ equation 
(JTJ) for any Borel set X results in 



J= J. (2) 

X J F{s,X) 

Equation (J2J) can be seen as follows: Let dD be oriented "anticlockwise" 
(direction of integration), then 



Since 



J = J+ J- J- J- 

<)D JX J l Pl {\r,r+s\) Jf(s,X) i 7p2 ([r,r+s]) 

(a) (b) 



r rP rP 

/ J = / J(%(s))ds= / J(v Jp{s) )ds 



f> , ff> 1 _ 

vj(jjE) Jp{s )ds = / -E(j,j) 7p{s) ds = 
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Fig. 1: Spacetime region D enclosed by integral curves and instanta- 
neous sets 



the contributions (a) and (b) to f gD J vanish: (a) = (b) = 0. 

For an interval ICE T integrals of the form (J2J) expressed in our chosen 
coordinates read as 



which is easily recognized as the standard quantum mechanical probability for 
"finding a particle" at time r in the spatial interval [mxn(x(X)), ma.x(x(X))}, 
as soon as p is integrable and normalized, i.e. f p(r, £)d£ = \I/*\I/(t, £)d£ 
1, which will further on be assumed to hold. Equation (J2J) can then be inter- 
preted as the conservation of probability along the flow lines of the Bohmian 
vector field. In other words, the amount of probability contained in X is 
the same as in F(s,X), for all s6K. Note that for a complete vector field 
v, we get a fibration of M. into the images of the integral curves or orbits 
of v . Equation (j3J) then induces a measure on the space of these orbits, 
independently of the choice of r, which will be argued below. 

Having this in mind, the amount of probability contained in the set of 
Bohmian orbits passing through quite arbitrary spacetime regions D, which 
need not be subsets of an instantaneous space, can be defined in a straight- 
forward and unambiguous manner: We simply take the probability Q of 
the intersection of all the flow lines, passing through D, with an arbitrary 
hypersurface E r . 




(3) 
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Let Ti be the projection n : M x M. — > M., (s,p) i— > p and let F T be the 
restriction of F onto ix E T . Then the mapping 



is the fibre projection of M. onto the instantaneous subspace S r along the 
Bohmian trajectories. If T T {D) is a Borel set in E T then we define the 
transition of the Bohmian vector field through D as 



JTt{d) 

With X T := xoT T delivering the x-coordinates of the fibre projection T T this 
yields 



Relation (J2J) secures that T(D) is independent of the choice of the instanta- 
neous subspace S T , or respectively of the choice of r. In what follows we will 
restrict ourselves to JF and X respectively. 

Remark 1 T does not induce a measure on the Borel subsets of A4, as 
T(D 1 UD 2 ) ^ T(D 1 )+T(D 2 ) in general, for DxnD 2 = 0. For regions D x C 
D 2 it is however guaranteed, that !F (Di) C JF (D 2 ); and therefore X (Di) C 
A'o(^2)> as a ^ ^ e trajectories passing through Di obviously also contribute 
to the transition through D 2 . I.e. for D\ C D 2 relation T(Di) < T(D 2 ) 
holds. 

3.3 Arrival time distribution 

Having the notion of transition through a spacetime region D, the definition 
of detection probabilities is obtained easily. For a spacetime region D C M. 
we define the subsets 



as indicated in figure 121 The transition T(D T ) of flow lines through D T then 
is the probability for a Bohmian particle, to have "arrived" in the region 
D before time r. If D is furthermore assumed to be the spacetime region 
occupied by a 100% efficient, purely passive detecting device, which "clicks" 
as soon as a particle, i.e. it's Bohmian trajectory, enters it, then the transition 
T(D T ) finds it's interpretation as detection probability up to time r. 




J. 




D T := {p E Dj t(p) < t} 



(4) 
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Fig. 2: Subsets D T of a spacetime region D 



Definition 2 For a spacetime region D with subsets D T Q) ; whose fibre 
projections Tq{D t ) are Borel sets for every r6l, the detection probability 



P D (r):=T(D T )= J J= J p(0,£)d£ 

is, as a function of time r G 1R, positive, monotonically increasing and 
bounded. 

The positivity is guaranteed because of the positivity of the probability 
density p = \E'*\1/, the monotonicity because of D T1 C D T2 for T\ < r 2 and 
remark [TJ P is bounded because the transition T(D) is bounded by 1 for 
any spacetime region D. 

Interpretatively the quantity N := lim -Pd(t) represents the overall prob- 

ability for a detection event in D. Now taking only that part of Bohmian 
world lines into account, which enter the spacetime region D at some time, 
and therefore produce a detection event, 4Pd(t) delivers the conditional 
probability for a detection event up to time r. 

Definition 3 Let Pp be the detection probability function of definition 
then by the function 

P C d--= y P p , T :R->[0,1] 
hm P D {r) 

the conditional arrival time distribution for a spacetime region D shall be 
denoted. 
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The lower index of P& and Pf, respectively will further on be omitted, if 
it is evident from the context, which spacetime region D is meant. 

4 Formulation without trajectories 

From now on we consider detectors occupying spacetime regions of the form 

D = {p <E M/x{p) <E [a, b], t(p) > 0} (5) 

which corresponds to a detector being at rest with respect to our inertial 
chart, occupying the interval [a, b], and being sensitive from time t — 
onward. For b — > a we get the point detector at x — a. 
Again we have the subsets 

D T = {pe D/ tip) < t}. 

Proposition 4 For the detection probability function P for a spacetime re- 
gion of type (0) the formula 

P( T ) = P(0) + max{/ a (s)/ < s < r} + max{-/ 6 (s)/ < s < r} (6) 

holds for times r > 0, with 

PS PS 

fa(s)= / j(t,a)dt and fb(s)= / j(t,b)dt 
Jo Jo 

being antiderivatives of the current densities at x = a and x = b, respectively. 



Proof: From remark^ we know that 

nmn(X (D T )) 

and also that 



P(r) = T(D T ) = J p(0, Ode > P(0) = T(D ) = J , 



mm(X()(D T )) < a and max(rYo(-D T )) > b 
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as Xq{Dq) = [a, b] C X Q (D T ), which allows the separation 

a b max(Ab(r> T )) 

p(r) = J p(o,e)d£+ J P (o,m+ J p(o,m 

mm(X (D T )) a b 

a mnx(X (D T )) 

J p(0,£)d£ + P(0)+ J p(0,£)d£. (7) 

min{X (D T )) b 

The non-crossing property of Bohmian trajectories secures that 
X ({(t,a)})<X ({(t,x)}), Vx>« 

and 

X ({(t,x)})<X ({(t,b)}), Vx<b. 
This allows the reformulation of (JJJ) into 

P(r) = P(0) 

+ f P(0,0d£ (8) 

Jmm(X (D«)) 

+ / P(0,0d£ (9) 

Jb 

with 
and 



-D" = {p G £> T / x(p) = a} 



D\ = {p e D T / x(p) = b} 



being the restrictions to the right, respectively left, edges of D T . Because 
of the positivity of the position density p, by applying Stoke's theorem to 
regions as illustrated in figure EH term (jSJ) can be expressed as 



Jmm(X (D2)) Ux Q ({(t,a)}) ) 

max < / j(s,a)ds > = max{/ a (£)/ < t < r} . 

K*<r [J J 



(10) 
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X ({(x,a)}) 

Fig. 3: Spacetime region 



Analogously term (jHJ) becomes 

nmax{X (D^.)) 
Jb 



p(0,Ode = max J / 
max < / j(s, b)ds 

0<t<r\J t 



*>({(*,&)}) 



= max{-/ 6 (t)/ 0<t<r}. (11) 

Equation (J7|) together with (JTUJ) and (fTTj) finally reads 

P(r) = P(0) + max{/ a (s)/ < s < r} + max{-/ 6 (s)/ < s < r} 
which completes the proof. ■ 

Taking the limit b — > a we immediately get the detection probability for 
point detectors: 

Corollary 5 For spacetime regions of the form 

D = {p G M/ x{p) = a, t(p) > 0} 
the detection probability function P can be expressed as 

P(t) = max{/ a (s)/ < s < r} + max{-/ a (s)/ < s < r} 

with s 

f a (s) = / j(t,a)dt 
Jo 

being an antiderivative of the current density at x = a. 
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The conditional arrival time distribution associated with the spacetime 
region (JSJ) 

pc ■= ~, — ^tt ■ R -> [0, 1] 
lim P(s) 1 J 

s^oo 

obeys P c (t) = for r < 0. P c (0) > holds for a < b, since p > 0. For 
r > P c is continuous and P c is non-decreasing. Thus P c is the distribution 
function of a Lebesgue-Stieltjes measure \i [llj. p can be separated into 
a point measure p s , which ascribes to every Borel set A C R the value 
P c (0) if G A, the value otherwise, and an absolutely continuous part 
/i ac , with respect to the Lebesgue measure. That is, there exists a Lebesgue- 
measurable density function 5 on R, which vanishes for r < and ascribes 
to every interval [a, b] the value 

ft ac ([a,b]) = P c (b) - P c (a) = I 5{s)ds 

J a 

if < a < b. fi finally reads as \i = fi s + fi ac 

The density 5 is needed for the calculation of the expectation values and 
variances of the arrival time, represented by the stochastic variable 

T A : R -> R, r (-»• r. 

If the relevant integrals exist, they can be obtained in the usual manner: 

poo /*oo 

(T A )» = (T A ), S + (T\ c = • P c (0) + / r • 5(r)dr = r ■ 5(r)dr 



and 



V,(T^) = ((T^) 2 ),-(T A )J. 



Remark 6 T7ie probability density 5 on ]0, oo[ takes the form 

8(t) = (lim Pit))'' ■ \j(r,a) ■ ( f a (r) - max {/„(*)} 

\t^OO J V 0<S<T 



J(r,6)-e -/ 6 (t)- max {-/ 6 (a)} 

0<S<T 



(12) 



•where B denotes the step function O(s) 



s < 

1 s > 



T7ie formulation for a point detector at x = a is simply achieved by re- 
placing b with a everywhere in Alty) . 
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This shows, that it is not enough to know the current density at a given 
instant r, rather the current density j{-,a) has to be known at all instants 
within the interval [0, r], in order to compute the probability density 5 of the 
arrival time distribution at a. Note however, that the Bohmian trajectories 
do not need to be known. The function 5 ■ lim P(t) is the cut off current 

t— »oo 

introduced in |B]. 



5 Examples 

The following numerical examples shall give an overview of the applicability 
of our treatment. 



5.1 Free evolution 

As a first example we choose a solution ^ to the free and parameter-reduced 
Schrodinger equation 

The parameter reduction is simply achieved by taking the t-coordinates in 
units of ^pr and the ^-coordinates in units of 1/q, where q is a characteris- 
tic wave number of the wavefunction (e.g. corresponding to a peak in the 
momentum space). We choose \I/ to be of the form 

%?(t,x) = )J^($(t,x;k ,x ) + $(t,x;k ,3xo) + $(t, x; k 0} 5x )) 

+ j| ( $ (^ x 'i ~ k o, -xo) + $(t, x; -ha, -3x ) + $(t, x; -k , -5ar )) 



with 



»(t, x; fr, xq) = ( d -) 7^ exp ( ^ ± ^ - *°>H 



and consequently initial data \&o = \l/(0, ■). This describes six differently 
weighted Gaussian wave packets moving with the same velocity in opposite 
directions. The evolution of this wave function for the chosen parameters 
ko — 5, d — 1 and x = —4 is illustrated in figure EJ 
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-40 -30 -20 -10 10 20 30 40 



Fig. 4: Free evolution of a one-particle wavefunction consisting of six 
Gaussians, moving to the right and left, respectively. 

The corresponding Bohmian trajectories, beeing solutions to the velocity 
vector field v, show the peculiar non- intersection property (figure (SJ). Even 
"free" particles change their direction of motion along their paths. The 
starting points of the worldlines are sampled according to \^o\ 2 - 




Fig. 5: Bohmian trajectories 
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We place a point detector D at x = —2.5. The corresponding arrival time 
distribution is gained with the method of proposition^] It shows areas of non- 
increasing arrival time probability, which correspond to times, during which 
already detected particles enter D for a second, third, etc. time, and therefore 
do not contribute to P(t) anymore. Figure EJ illustrates this phenomenon and 
the technical procedure of proposition 
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Fig. 6: f a (dotted line), max{/ Q (s)/ < s < t}, -max{-/ a (s)/ < 
s < t} (upper and lower dashed lines), the detection probability 
function P (dashed-dotted line) and the conditional arrival time 
distribution P c (solid line). 



5.2 The potential barrier 

Terms as arrival-, delay- or dwell-times, etc in the literature are often con- 
nected to scattering states of one-dimensional potential barriers. For reviews 
of the subject see [TJJUHHIll . We apply equation © to a situation like that. 

Now the evolution of our wavefunction \l/ is given by the Schrodinger 
equation 

id t m = ~dl + -(e(x-a)-e(x-b)) ^ 

with a < b and G denoting the step function. The parameter reduction in 
this case can be achieved by taking the t-coordinates in units of t^t and the 
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x-coordinates in units of 7^=^ , where Vo is assumed to be the height of the 
potential barrier. 



0.2 




Fig. 7: Evolution of a Gaussian wave packet impinging on a potential 
barrier. 

The initial gaussian wave packet \l/o is placed sufficiently far to the left 
of the potential barrier (for the purpose that at time t = the interference 
effects due to the barrier are negligible). \I/ is moving towards the barrier. 
The evolution of the wavefunction is illustrated in figure Figure |H] shows 
the corresponding Bohmian trajectories sampled according to l^ol 2 - 

We consider an extended detector occupying the spacetime region of the 
barrier from time t = onwards. As the wave packet is placed sufficiently 
far to the left of the barrier, the probability of a detection event at t = is 
negligible. Figure |H1 shows, that is positive and monotonically increasing, 
which corresponds to a positive current density at the right edge of the 
potential barrier. Therefore /& does not contribute to an increase of the 
detection probability. As trajectories enter the detector from the left, f a at 
the left edge of the barrier increases up to a time r r when the first trajectory 
is reflected before entering the detection region. From r r onwards then part of 
the trajectories inside the barrier return to be finally reflected and therefore 
produce a negative current density at x — a, which leads to a small decrease 
of f a . The resulting detection probability function P and the conditional 
arrival time distribution P c are indicated in figure M by the dashed-dotted 
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-60 -40 -20 20 40 60 

Fig. 8: Bohmian trajectories 
and solid lines, respectively. 




10 20 T 30 40 50 60 

r 



Fig. 9: f a and fa (upper and lower dotted lines), the detection proba- 
bility function P (dashed-dotted line) and the arrival time dis- 
tribution P c (solid line) for the Gaussian wave packet impinging 
on a potential barrier. 
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5.3 The potential step 



As a third numerical example we take the case of total reflection at a potential 
step at x = 0. The wavefunction \l/ is now the solution of the Schrodinger 
equation 



--dl + -Q(x) 
2 x 2 K ' 



with again denoting the step function. The parameter reduction can be 
achieved analogously to the case of the potential barrier. Again the initial 
data \I/o of \I/ is taken to be a Gaussian wave packet placed to the left of the 
potential step at t — (figure ITUj) . 




Fig. 10: Evolution of a Gaussian wave packet being totally reflected at 
a potential step at x = 0. 



This time we consider the situation, that at the time the detector is 
activated, some of the trajectories are already located inside the detection 
region D, located between x = a and x = b in front of the potential step, as 
illustrated in figure [TT] 

The detection probability at time t = thus takes a value noticeably 
different from 0. As we have the case of total reflection, all the trajectories, 
which initially started to the right of b, eventually turn back and pass the 
detector at a later instance. Therefore almost all the trajectories pass D at 
some time and the limit lim P(s) is approximately 1. The detection prob- 
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-50 -40 -30 -20 -10 10 

a b 

Fig. 11: Bohmian trajectories 

ability function and the conditional arrival time distribution thus coincide 
and are indicated by the solid line in figure IT21 
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-0.6 - 
-0.8 - 

-1 - 

10 20 30 40 50 60 70 80 90 100 

Fig. 12: f a and /& (dotted lines), max{/„(s)/ < s < t} and 
— max{— fb(s)/ < s < t} (upper and lower dashed lines), 
the detection probability function P and the conditional arrival 
time distribution P c (solid line). 
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